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Abstract 

We report a determination of the antiproton magnetic moment, measured in a three-body system, independent of 
previous experiments. We present results from a systematic study of the hyperfine (HF) structure of antiprotonic 
helium where we have achieved a precision more than a factor of 10 better than our first measurement. A comparison 
between the experimental results and three-body quantum electrodynamic (QED) calculations leads to a new value for 
the antiproton magnetic moment p v s = -2.7862(83)//ai, which agrees with the magnetic moment of the proton within 
2.9 x 1(T 3 . 
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Antiprotonic helium (pHe + ) is a three body system 
consisting of an antiproton, a helium nucleus and an 
electron (p, He ++ , e~) HI Etl- It is formed by stop- 
ping antiprotons in a helium medium. Because of its 
mass, the p occupies states of principle quantum num- 
ber n ~ 38 with the highest probability, while the e~ re- 
mains in the ground state. The vast proportion of newly 
formed pHe + atoms proceed rapidly (within nanosec- 
onds) to an ionised state by Auger excitation of the elec- 
tron. The antiproton then annihilates almost instanta- 
neously (within picoseconds) with one of the nucleons 
in the helium nucleus due to the overlap of their wave- 
functions. However, a small proportion (~ 3%) occupy 
circular states n ~ /, where I is the total angular momen- 
tum quantum number. Auger decay is then supressed by 
the large ionisation energy (~ 25 eV) and degeneracy 
is lifted due to the presence of the electron |4]. These 
states become relatively long lived (metastable) because 
the only decay channel available to them is the radiative 
one, and they cascade from («, I) — > (n — 1, 1 — 1) with 
typical energy level spacings ~ 2 eV and lifetimes of 
~ 1.5 ps. 

Shortly after its discovery, in 1991 J5h, laser spec- 
troscopy measurements were performed on various lev- 
els of the pHe + cascade 0. S Compared with 
three-body Quantum Electrodynamic (QED) calcula- 
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tions these measurements were used to determine the 
antiproton-to-electron mass ratio (mp/m e ). Over the 
years, improvements to the experimental system have 
increased the precision to 2 ppb fioll . which is one of 
the best tests of CPT invariance in the baryon sector. In 
1997 the hyperfine structure was first revealed by laser 
spectroscopy 111 111 and, in 2002, measured via a laser- 
microwave-laser technique 1 1 211 . 

This paper concludes a systematic study, commenced 
in 2006 |13fl , where the error has been reduced by more 
than a factor of 10 over the first measurement 11211 and a 
new value for the antiproton spin magnetic moment p p s 
has been determined by comparison with QED calcula- 
tions. 

The hyperfine (HF) splitting |Ql of pHe + arises from 
the coupling of the e~ spin § e with the p orbital angular 
momentum \L\ ~ 35h and results in a doublet structure 
of the order vhf = 10 - 15 GHz. The interaction between 
the antiproton spin S"p with S ' e and £ causes a superhy- 
perfine (SHF) splitting of size v SHF = 150 - 300 MHz. 
A schematic of the energy level structure is presented in 
Fig. El. 

The theoretical framework for the level splitting has 
been developed by two separate groups IhI [l5i [la, 17 1 



which all use the same Hamiltonian, first derived by 
Bakalov and Korobov II 1411 . but different variational 
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Figure 1: (a) Schematic view of the level splitting of 
p~He + for the (n, I) — > (n — 1, I + 1) electric dipole 
transitions. The laser transitions / + and f~, from the 
parent to daughter states, are indicated by straight lines 
and the microwave transitions, between the quadru- 
plets of the parent, by wavy ones. For this experiment 
(n,L) = (37, 35) and (n',L r ) = (38, 34). (b) Laser reso- 
nance profile demonstrating the two sharp peaks and HF 
laser splitting Af = f~ - f + ■ Although there are four 
SHF laser transitions only the HF ones can be resolved 
in this experiment. 



methods to extract the energy eigen values. 

The HF doublet is characterised by the quantum num- 
ber F — £ + S e and the SHF quadruplet by /= F + SW = 
L + § e + 3p. An electron spin flip transition can be in- 
duced by an oscillating magnetic field, resulting in two 
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Bakalov and Widmann II18H indicate the sensitivity of 



certain states on /j p s and the precision required to im 



prove its value over the most precise measurement 11911 . 
Some of these states are not practical due to limitations 
in laser capability. Others, like the («, L) = (39, 35) 
state, which are within laser capabilities, have a HF laser 
splitting Af = 0.5 GHz, of the same order as the Doppler 
broadening A/b = 420 MHz. The previously measured 
(37, 35) state remains the best candidate for a precision 
study because there is an easily stimulated laser transi- 
tion between the (n, L) = (37, 35) and (38, 34) states 
with Af = 1.75 GHz. It is therefore possible to indi- 
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Figure 2: (a)Schematic to show timing between the mi- 
crowave envelope, the two laser pulses and the resulting 
annihilation peaks. The microwave envelope starts be- 
fore the first laser and ends after the second. However 
before firing the first laser there is no population asym- 
metry and thus no observable transition can occur. Af- 
ter firing the second laser no further measurements are 
made so the effective microwave pulse length is equal 
to the delay between the two lasers T. (b) Schematic of 
the target region. 



vidually resolve the F ± states, demonstrated in Fig.QJ), 
since Afo = 320 MHz. The parent state is also relatively 
highly populated, containing some 0.3% of the antipro- 
tons stopped in the target 12011 . 

The laser spectroscopy experiments employ a tech- 
nique by which the annihilation decay products are 
detected J2|]. A sharp prompt peak is first observed, 
where the majority of states annihilate within picosec- 
onds of formation, then an exponential tail, where the 
metastable states cascade more slowly towards the nu- 
cleus. This tail constitutes the background. 

A narrow band pulsed laser is scanned over the re- 
gion of an expected transition between a radiative decay 
dominated parent state and an Auger decay dominated 
daughter. Because the daughter state is relatively short 
lived (~ 10 ns), resonance is indicated by a sharp peak 
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Figure 3: Microwave frequency profiles averaged from 
scans at a common pressure p — 150 mbar and laser de- 
lays T = 200 ns, 350 ns and 500 ns. The broadening 
is due to the Fourier transform of the rectangular mi- 
crowave pulse of length T. 



against the background at the time of the pulse. The ra- 
tio of the peak area to this background (peak-to-total) 
indicates the size of the population transferred. 

The HF splitting measurement method is illustrated 
in Fig. [2a] where the first laser pulse remained fixed 
to the / + transition between the (37, 35) and (38, 34) 
states, creating a population asymmetry. A microwave 
pulse, if on resonance with either y+p or v HF , then trans- 
ferred the population from J~ + or J to refill the F + 
state. A second laser pulse was tuned to the same / + 
transition and fired with a delay T = 200 - 500 ns from 
the first, which measured the population transfer. Plot- 
ting the peak-to-total of the second laser induced annihi- 
lation peak as a function of microwave frequency yields 
the two HF transitions as distinct peaks 113 1. 

The experiment was carried out at CERN's Antipro- 
ton Decelerator (AD) which provided a pulsed beam 
of(l-3)xl0 7 antiprotons with length ~ 200 ns and 
energy 5.3 MeV. Every 90 - 120 s, such a pulse was 
stopped in a helium gas target at a temperature of 6.1 K 
and pressures p = 150 - 500 mbar (number density 1.7- 



6.2 x 10 cm ). The antiproton annihilation products 
passed through one of two Lucite plates either side of 
the target where their Cherenkov photons were detected 
by photomultipliers (PMT) 12 111 and the signal was dis- 
played on a digital oscilloscope. The PMTs were gated 
off for the p arrival so that only the 3% metastable tail 
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Figure 4: The difference A exp -th = v exp - v tn for vSgJbe 



tween experiment and the closest theory = MHz II 1511 . 
The second theory lfl7ll is another 300 kHz less. The 
experimental value of v^p is shown as a solid triangle 
(A) and v HF as an empty tr iang le (A). The estimated 
theoretical error is 1.3 MHz ill 80 and therefore too large 
to be shown on the scale of these graphs, (a) Pressure 
dependence, where the point at 250 mbar is the average 
of two power dependent measurements from [13]. The 
points represented by the solid circle (•) and empty 
circle (O) shown at p — mbar are v^p and Vp, F re- 
spectively, (b) Power dependence measured at constant 
pressure p — 150 mbar the average of which constitutes 
the point at p = 150 mbar in (a). 



was recorded. 

The pulse-amplified continuous wave (CW) laser sys- 
tem | lol 22 ] was constructed by splitting a CW laser of 
wavelength A ~ 726.1 nm into two seed beams. These 
were each pulse amplified by a NG:Yag laser and three 
Bethune dye cells, the second delayed by a time T af- 
ter the first. The pump beams were stretched so that the 
pulse lengths of the two lasers were 18 ns and 13 ns lfl3tl . 



Table 1 : A list of the experimental data including parameters pressure p, laser delay T, microwave power P, number 
of shots from the AD p shots, HF transition frequencies v^ F , and peak widths r ± . 



p (mbar) 


T(ns) 


P(W) 


p shots 


v+ F (GHz) 


T + (MHz) 


v HF (GHz) 


T_ (MHz) 


150 


200 


15 


1070 


12.896 45(12) 


3.84(8) 


12.924 30(14) 


4.70(3) 


150 


350 


5 


1028 


12.896 709(78) 


2.76(8) 


12.924 579(59) 


2.21(4) 


150 


500 


3 


2236 


12.896 688(48) 


1.68(6) 


12.924 470(41) 


1.48(4) 


250 [13J 


350 


5 


2938 


12.896 651(35) 


2.24(2) 


12.924 431(35) 


2.41(3) 


250 [13J 


500 


3 


230 


12.896 53(12) 


2.55(7) 


12.924 446(65) 


1.65(5) 


500 


350 


5 


1844 


12.896 525(80) 


2.06(5) 


12.924 446(99) 


2.01(7) 



slightly longer than the Auger lifetime of the daughter 
state, to achieve a high depopulation efficiency. 

The microwave apparatus was similar to that de- 
scribed in Sakaguchi et al. 12311 and a schematic is dis- 
played in Fig. [2b] The microwave pulse was synthe- 
sised by a vector network analyser (Anritsu 37225B) 
referenced to a 10 MHz satellite signal (HP 58503B) 
and amplified by a travelling wave tube amplifier (TMD 
PTC6358). A waveguide carried the pulse to a custom 
made stainless steel cylindrical cavity, with central fre- 
quency ~ 12.91 GHz, which provided the desired shape 
for the field (TMno mode) at the target. Steel meshes 
(92% transparency) covered both ends of the cavity so 
that antiprotions and the two laser beams could enter the 
target from opposite directions. The cavity was over- 
coupled to the waveguide to achive a broad frequency 
range ~ 100 MHz. A mu-metal shell surrounded the 
target region to protect from external magnetic fields. 
Indeed the field measured in three dimensions within 
the target was B < 0.03 G. 

Previously, different choke positions of a triple-stub- 
tuner were used to match the impedance of the waveg- 
uide to that of the cavity for a range of frequencies 112311 . 
This time, a constant microwave power P was produced 
at the target by firing a predetermined signal strength 
down the unmatched waveguide. Most of the signal was 
reflected and dumped to a 50 O terminator by a three- 
way circulator. This removed standing waves from 
the system and allowed the relatively small amount of 
power absorbed by the cavity to be controlled to within 
1 dB over the frequency range. The power was mon- 
itored by an undercoupled pickup antenna situated op- 
posite the waveguide. 

Table Q] shows a summary of all data measured for 
this experiment. The line shape was determined in Pask 
et al. II 1 3fl to be 
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Data measured with the same p and T were fitted simul- 
taneously with common parameters for height, width 
and central frequency. Two data sets were systemati- 
cally examined: 1) Microwave power dependence, and 
2) Pressure dependence 

1) The ac Stark effect shifts the El transitions by less 
than one part in 10 9 ifioll . Its equivalent, the ac Zee- 
man shift of the Ml transitions is far weaker and there- 
fore far too small to be resolved. A power dependence 
measurement was nevertheless examined for a complete 
understanding of the systematics. At a constant pres- 
sure p = 150 mbar, resonance profiles were measured 
with various laser delays T = 200 ns, 350 ns and 500 ns 
at microwave powers carefully chosen to achieve a n- 
pulse lEl]. P = 15 W, 5 W and 3 W, respectively. For 
illustrative purposes the average of these scans is shown 
in Fig. [3] where the the dominating broadening effect 
was due to the Fourier transform of the rectangular mi- 
crowave pulse of length T [25]. Figurel4blshows that no 
power dependent trend was observed, therefore all data 
measured at common target densities could be averaged. 

2) In 2006 Korenman predicted a collisional shift of 
A = dv/dp = 0.3 kHz/mbar lE6Tl and, for the first time, 
a resolution has been achieved to examine this. High 
statistic microwave resonance profiles were measured 
at p = 150 mbar and 500 mbar. Previous measurements 
had been made at p = 250 mbar 111 311 . 

The results are displayed in Fig. [4a] and clearly 
show that, if any such density dependence existed, 
the trend may have the same magnitude but op- 
posite sign. Extrapolating the points to zero den- 
sity can be performed with lines of average gradient 
A = -0.26 ± 0.2 kHz/mbar. Such a large error neither 
confirms nor precludes conclusively the existence of a 
shift. The two fits in Fig. [4a|are shown with the same 
gradient because no density shift was observed for the 
difference between the transitions Avhf : 
Fig. [5] in accordance with predictions. 

A more recent calculation 12711 . based on mea- 



'HF 



surements of the collisional broadening 12511 . predicts 
A = -0.048 kHz/mbar which is less than the experi- 
mental precision. As a negative slope contradicts the 
predictions, the average of each measurement was taken 
but the error used was that of the extrapolation to zero 
density from the fit. 
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Figure 5: The difference for Avhf between experiment 



and theory II15L Il7ll as a function of target pressure. 
The experimental values are shown as squares (■). The 
point shown at x — and represented by a circle (•) 
represents the average of the total data. 

Other systematic effects that influence the measure- 
ment include external magnetic fields, precision of the 
microwave frequency source, shot-to-shot microwave 
power fluctuations and variances in the laser position 
and fiuence from day to day. However these effects 
have been determined to be far smaller than the shot- 
to-shot fluctuations of the antiproton beam. Data was 
measured over a long period to reduce these drift effects 
and variations in the p intensity have been reduced by 
normalising the second laser induced annihilation peak 
with the first (proportional to the number of antiprotons 
captured). Despite these considerations the reduced chi- 
squared Xied of the fit was Xmd ~ 3. To adjust for this 
the error bars were inflated by VVred ~ 1.7. 

Bakalov calculated that a broadening of v* F due 
to an external magnetic field occurs at a rate of 
T ± -5.6 MHz/G H. The similarity between the 
Fourier transform of the microwave pulse and the spec- 
tral line widths [25] confirms that the target region was 
well shielded during the experiment. Due to referenc- 
ing to a 10 MHz GPS receiver, the precision of the fre- 
quency source is several orders of magnitude less than 
the resolution of this experiment. Thefore the statistical 
errors are much greater than the systematic. 

The individual transition frequencies have a negligi- 



ble dependence on jj^. However Avhf is directly propor- 
tional to this value. The predicted density shift for Avhf 
is far smaller, A = 0.003 Hz/mbar l27ll than the precision 
of this experiment. If this is the case, the total split- 
ting can be calculated from the difference between each 
pair of transitions that were measured at common den- 
sities Avhf = ( v hf _ v hf )/^> ratner than the differ- 
ence between the sum of each transition measurement 

Avhf = (Zf v hf _ 2f v hf -)/^> wnere ! is tne index of 
a measurement and N - 3is the total number of density 
dependent measurements. Figure [5] displays Avhf as a 
function of target pressure compared to the two most 
recent theories. 

Fitting a first order polynomial, results in a gradient 
almost half that of its associated error, A = 0.24 + 0.37 
kHz/mbar, so the above holds and the data can be aver- 
aged to obtain a final value of Avhf- Table|2]presents the 
data for the recent and previous experiments compared 
to the two most up to date theories. 

This work demonstrates the completion of a sys- 
tematic experimental study on the HF splitting of the 
(37, 35) state of pHe + . The experimental error of v* F 
has been reduced by a factor 20 less than that of the 
theoretical calculations and, although A exp _ t h = 300 - 
600 kHz, it is well within the estimated theoretical er- 
ror 1.3 MHz. The experimental precision for Avhf has 
reached that of theory and has been improved by more 
than a factor of 10 over the first measurement 111 211 . 
There is a two sigma agreement between theory and ex- 
periment. 

The sensitivity S of Avhf on /j^ for the (37, 35) state 
is S = dE/d[£ = 10.1 MHz/jUiv 03, where fx N is the 
nuclear magneton. Thus the magnetic moment can be 
determined to be: 



-2.7862(83)^ 



(1) 



where the uncertainty has been calculated by adding 
A exp -th with the errors of theory and experiment in 
quadrature, resulting in a one sigma error, slightly less 
than the value determined by Kreissl et al. 11911, while 
the deviation from the magnitude of the proton spin 
magnetic moment, yuf = 2.792847351(28), is similar 
but opposite in sign. 

The absolute values for the magnetic moments of the 
proton and antiproton are in agreement within 



^ l / sl = (2.4 ± 2.9) x 10- 3 . 



(2) 



The limit of experimental precision has been reached 
for the (37, 35) state. The study of other states with 



Table 2: Experimental data compared with three-body QED predictions, where v^ F are the HF transition frequencies 
and Avhf is the difference between v HF and v^ F . The quoted theoretical errors have been estimated by Bakalov and 



Widmann (1811 





v^ F (GHz) 


v HF (GHz) 


Av HF (MHz) 


This work 
2002 [12] 


12.896 641(63) 
12.895 96(34) 


12.924 461(63) 
12.924 67(29) 


27.825(33) 
28.71(44) 


Korobov [15] 
Kino [17] 


12.896 3(13) 
12.896 0(13) 


12.924 2(13) 
12.923 9(13) 


27.896(33) 
27.889(33) 



larger Avhf could potentially increase the precision 
but the system cannot generally be improved due 
mainly to uncontrollable fluctuations of the AD beam. 
Preparations are underway to measure the HF splitting 
of p 3 He + which, because of the additional helion spin, 
provides a more thorough test of the theory but yields 
no further information on the magnetic moment. The 
theorists are currently working on an a 6 calculation but 
a significant change in Avhf is not expected. 
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